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Abstract 

Density matrix for N-qubit symmetric state or spin-j state (j = N/2) is expressed in terms of the 
well known Fano statistical tensor parameters. Employing the multiaxial representation [l], wherein 
a spin-j density matrix is shown to be characterized by j(2j+l) axes and 2j real scalars, we enumerate 
the number of invariants constructed out of these axes and scalars. These invariants are explicitly 
calculated in the particular case of pure as well as mixed spin-1 state. 
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1 Introduction 

The problem of enumeration of local invariants of quantum state described by a density matrix p is 
important in the context of quantum entanglement. Nonlocal correlations in quantum systems reflect 
entanglement between its parts. Genuine non- local properties should be described in a form invariant 
under local unitary operations. Two N-qubit states are said to be locally equivalent if one can be trans- 
formed into the other by local operations, i.e., p' = UpW where U G SU(2) xN and the two quantum 
states p and p' are said to be equally entangled. 

A general prescription to identify the invariants associated with a multiparticle system have been 
outlined by Linden et al., \2\- Well known algebraic methods for generating invariants already exists 
in literature (3||6j. Williamson et al., [7] have presented a geometric approach for constructing SU(2) 
and SL(2,C) invariants. Makhlin [8] has presented a complete set of 18 local polynomial invariants of 
two qubit mixed states and demonstrated the usefulness of these invariants to study entanglement. As 
the number of subsystems increases, the problem of identifying and interpreting independent invariants 
rapidly becomes very complicated. Usha et al., |9j have shown that a set of 6 invariants which is a 
subset of a more general set of 18 invariants proposed by Makhlin [8] is sufficient to characterize the non 
local properties of a symmetric two qubit system. We also focus on symmetric two qubit states as the 
problem of identifying independent invariants become easier. Our approach makes use of the geometrical 
multiaxial representation of an arbitrary spin-j density matrix [I] which is completely characterized by a 
set of j(2j+l) axes and 2j real positive scalars. 
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The paper is organized as follows: In Sec. 2, the decomposition of a density matrix in terms of 
well known Fano statistical tensor parameters is presented. Discussion of multiaxial description of the 
density matrix using Wigner-D matrices is outlined and the invariants associated with N-qubit symmetric 
state are enumerated in section 3. In Sec. 4, We have explicitly calculated the invariants of two qubit 
symmetric mixed as well as the most general pure state. To make our task easier, we have considered 
the Special Lakin Frame which is widely used in nuclear reactions. 



2 Symmetric subspace 

Here we are interested in the set of N-particle pure states that remain unchanged by permutations of 
individual particles. Symmetric states offer elegant mathematical analysis as the dimension of the Hilbert 
space reduces drastically from 2 N to (N + 1), when N qubits respect exchange symmetry. Such a Hilbert 
space is considered to be spanned by the eigen states {\j,m);—j < m < +j} of angular momentum 
operators J 2 and J z , where j = y. Analyzing a general state of N-particle spin-1/2 system represented 
by the density matrix of dimension 2 N x 2^ is difficult because the system's Hilbert space increases 
exponentially with the number of qubits N. Fortunately, a large number of experimentally relevant 
states possesses symmetry under particle exchange and this property allows us to significantly reduce the 
computational complexity. Completely symmetric systems are experimentally interesting, largely because 
it is often easier to nonselectively address an entire ensemble of particles rather than individually address 
each member and it is possible to express the dynamics of these systems using only symmetry preserving 
operators. The symmetric subspace therefore provides a convenient, computationally accessible class of 
spin states. Specifically, if we have N two level atoms, each atom may be represented as a spin- 1/2 
system and theoretical analysis can be carried out in terms of collective spin operator J = TjY,^ =l <T a . 
Here a a denote the Pauli spin operator of the ath qubit. The standard expression for the most general 
spin-j density matrix in terms of Fano statistical tensor parameters t q ' s is given by 

k — q — — k 

where r q (with Tq = I ,the identity operator) are irreducible tensor operators of rank 'k' in the 2j+l 
dimensional spin space with projection 'q' along the axis of quantization in the real 3-dimensional space. 
The T q satisfy the orthogonality relations 

Tr(rfr^) = (2j + l)S kk rS qq/ . (2) 



Here the normalization has been chosen so as to be in agreement with Madison convention 10 . The 



spherical tensor parameters t q 's which characterize the given system are the average expectation values 
given by t q = Tr j^ g - . Since p is Hermitian and = (— l) q T k q , t q satisfy the condition 

# = (-!)*** (3) 
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The spherical tensor parameters t q s have simple transformation properties under co-ordinate rotation 
in the 3-dimensional space. In the rotated frame t q 's are given by 

(t k q ) R = £ Dj s foM)# , (4) 

q = — k 
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where D*, (<fi,8,ifi) denote Wigner-D rotation matrix described by Euler angles ((f), 9, ip). 

3 Multiaxial description of density matrix 

It has already been shown [T] that a spin-j density matrix is characterized by j(2j+l) axes and 2j real 
positive scalars. For the sake of completeness, we reproduce it here; In general t±, can be made zero for 
any k by suitable rotation, i.e., 

+k 

(t k ±k ) R = = £ D k qfj±k (ct>,e^)t k ql . (5) 

q'=-k 

Using the well known Wigner expression for the rotation matrix D k , the above equation can be 
written as 

2k 

(t k ±k ) R = = [±% (9/2)} 2k exp[i(0 + VO] £ ^ » (6) 

r=0 

where the complex variable Z = cot(9 /2)e~ 1 ^ in the case of = and Z = tan(9 /2)e -4 ^ +7r ) in the 

case of = 0. The expansion coefficients C r in the polynomial are the same in both the cases and is 

i 

given by C r = = *r-fe • By solving the above polynomial equation, one can get 




in general two sets of k-coordinate frames in which (t± k ) = 0. Explicitly, if t k = in a coordinate system 
where Z axis is directed along (9, <j)) in the laboratory, t k _ k = in a coordinate system where Z axis is 
directed along (it — 9, (j) + tt). One set is obtained by the other by inverting the Z-axis. Therefore it is 
sufficient to enumerate the k independent solutions Qi(9i, 4>i), i=l,2,....k which constitute any arbitrary 
t k as a spherical tensor product of the form 

t k = r fc (...((Qi 8) Q 2 ) 2 ® Q 3 ) 3 ® ••••) fc ~ 1 ® Qkf q , (7) 

where 

(Qi®Q2)l = ^C(Uk;q 1 q2q)(Q 1 ) <ll (Q2) g2 I (Q)q = yyl r i g (M)- (8) 

Here C(11A;; qiq2Q) is the Clebsch Gordan Co-efficient and <^>) are the well known spherical harmon- 

ics. If one of the Qi 's is inverted, sign of equation(7) is changed. Hence it is possible to choose k axes 
Qi's , i=l,2,....k in such a way that r k is always positive. Each axis requires two independent parame- 
ters (9, 4>) to characterize it, hence the k axes together with the overall multiplicative factor account for 
exactly (2k+l) real parameters needed to specify a spherical tensor t k satisfying equation (3). Thus any 
spherical tensor of rank k can be represented geometrically by a set of k vectors Qi on the surface of a 
sphere of radius r. Consequently, the state of a spin-j assembly can be represented geometrically by a set 
of 2j spheres, one corresponding to each value of k, k=1....2j, the kth sphere having k vectors specified 
on its surface. 

Since (Qi(9i,4>i)(g)Qj(9j,(f)j))Q is an invariant (i 7^ j), one can construct in general ( ^ ^ ^~ ^ \ in- 
variants from j(2j+l) axes. Together with 2j real positive scalars, there are [ ^ ^ ' ^ J + 2j invariants 
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characterizing spin-j density matrix. Thus using this multiaxial parametrization of density matrix, we 
enumerate the total number of SU(2) invariants characterizing a spin-j density matrix . Let us consider 
the example of two qubit symmetric state for a detailed discussion. 



4 Invariants of two qubit symmetric state or spin-1 state 
4.1 Pure spin-1 state 



Consider the direct product \ip\) <8> \tp2) of two spinors in the qubit basis as 



1^12) 



cos- 



cos- 



( cos^cos^ 
cos^sin^e^ 2 



\ 



sin^-sm^-e^ 1 ^ 



2 



(9) 



< 0i,2 < 7T, < 0i,2 < 2ir . In the symmetric angular momentum subspace |11), |10), |1 — 1), the 
combined state will have the form 



1^12) 



sym 



( 

V 



02 



COS -j COS 2 

^(cos^sin^e^ 2 + sin^cos^e^ 1 ] 
sin^sin^e^ 1 ^ 



\ 
J 



(10) 



Since the two directions (0i,0i), (02,02) associated with the above two spinors define a plane, we 
choose this to be the xz-plane with respect to a frame xqUqZq with zq being the bisector of the above 
two directions. Thus the azimuths of the above two directions (0i,</>i), (02,02) with respect to xq are 
respectively and n. If the angular separation between the two directions is 20, then the state has 
the explicit form 



= \cos -111) $ — sin -1 
Vl + cos 2 9 l 2 1 ,z ° 2 1 



1) 



(11) 



The density matrix corresponding to the above state is given by 

/ 



(1 + cosH 



cos% 




-sin 2 § cos 2 \ \ 




\ -sin 2 \cos 2 \ 



sin 



4(9 



(12) 



/ 



Comparing equation (12) with the standard representation of the density matrix 



Tr(p) 



1 + 



3+1 + 3. 
2 J 0+ ^2 



V 



Vst 2 



1-V2t 2 



3 (t\ ~ t\) 



\ 



V3t 2 _ 2 



3 A + io_ , 
2 r + ^/ 



(13) 



we get the non-zero t^'s to be 



\/6cos0 
1 + cos 2 9 



t 2 



1 2 2 \[?>sin 2 d 
, t 2 = £_ 2 = — - cos20 ^ . 
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Since t\_ x = , zq itself is the axis (Qi) associated with t . As tj = ri(Qi)o , 

n = S- . (14) 
(Qx)l 

Solving for the polynomial equation (6) for t 2 , we get (Q 2 )l = \f^Yq(9,0) and {Qz)\ = \pfYq(6, it). 
Hence 

t 2 t 2 

r 2 = . h . = . h , . (15) 

(Q 2 ® (Q 2 ® Q 3 )| 

The invariants associated with the most general pure spin-1 state are 

h = n , h = r 2 , h = (Qi ® 2 )o ,i4 = (Qi ® Qz)l , ^5 = (Q 2 ® Qs)o ■ (16) 

Explicitly, 



cosv 

-*1 = T~l 97i > -'2 = — ^ , J3 = ^4 = 7^ , ^5 



cos2<9 



l + cos 20>- i + cos 2^'- - V3 V3 ' (17) 

It is clear from equation (11) that the state \tp) is separable for 9 = and it . Hence the invariants in the 
case of pure spin-1 separable states are 

ft t ^ t i 1 t 1 

Jl -V2 ' /2 -^' /3 - l4 -^' /5 --^- 

4.2 Mixed spin-1 state 

Consider the example of channel spin-1 system which plays an important role in nuclear physics 
experiments like hadron scattering and reaction processess [IT- 15 . A beam of nucleons colliding with a 
proton target provides such an example. If both the beam and the target are prepared to be in mixed 
states, then the corresponding density matrices are given by 

p(i) = ^[I + a(i)-p(i)] = lY, t «( i ) T ?^ i = 1 » 2 - ( 18 ) 

k,q 

where p(i) are the polarization vectors and a(i)'s are the Pauli spin matrices. 

The combined density matrix is the direct product of the individual density matrices 



p c = p{l) ® p(2) . (19) 

Eventhough the combined density matrix is a direct product of individual matrices, in this case, 
entanglement appears due to the projection of the combined density matrix onto the desired spin-1 
space. While solving this problem the Special Lakin frame (SLF) which is widely used in studying 
nuclear reactions is considered : Choose zq to be along p(l) + p(2). Since p(l), p(2) together define a 
plane in any general situation, we choose xo to be in this plane such that the azimuths of p(l) , p(2) 
with respect to xq are respectively and it. The yo axis is then chosen to be along zq x xq. The frame so 
chosen is indeed the special Lakin frame (SLF) as here i±i = and t 2 = t 2 _ 2 ■ Choose a simple case of 
|p(l)| = \p(2)\ = p, then we get i 2 ^ = in SLF. The density matrix so obtained for spin-1 mixed system 
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Figure 1: xoUqZq frame with mean spin direction zq as the bisector of two directions p(l) and p{2) 



in the symmetric subspace |11), 1 10) and |1 — 1) is 



/(l +pcos6f 



(3 + p 2 cos29) 



\ —p 2 sin 2c 



—p 2 sin 2 8 \ 



1-p 2 



(l-pcos9) 2 J 



(20) 



Observe that when p=l, the mixed state density matrix is exactly the same as that of pure state density 
matrix as given by equation (12). Comparing the above density matrix with the standard form (equation 
(13)), we get the non-zero vl's as 



2\^6pcos9 
(3 + p 2 cos29) 



t 2 



V2p 2 (l + cos 2 6) 
(3 + p 2 cos26) 



i L 2 



vV'" ; ' 



sin 



(3 + p 2 cos29) ' 



Solving for the polynomial equation (6) for t , t 2 , we obtain Q± = zo , Q2 = p(l) and Q3 = p(2) . Thus 
the invariants associated with the most general mixed spin-1 state are found to be 



2V6p\cos0\ 
1 ~ (3 + p 2 cos29) ' 2 



2^3p 2 



(3+p 2 cos26) 



, I3 = I4 = cosO , h 



cos20 



(21) 



Note that in both pure as well as mixed state, ^3 = I4 



cosd 
" V3 ' 



h 



cos26 
' V3 



For p = 1 and 9 = 0, it , the 



state is separable as in the case of pure state. For p < 1, the state is seperable for a range of values of 9. 
It is observed that as p decreases, the region of 9 for which entanglement appears also decreases [16]. 
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5 Conclusion 

We have considered symmetric N-qubit density matrix expressed in terms of Fano statistical tensor 
parameters. Making use of the well known multiaxial decomposition of the density matrix we have 
enumerated SU(2) invariants of the most general symmetric state. Considering the special case of two 
qubit symmetric state we have explicitly computed 5 invariants which form a complete set. Our framework 
can be applied to enumerate a complete set of invariants of any qudit state. The study of the relationship 
between various measures of entanglement and our complete set of invariants is in progress. 
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